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ABSTRACT 



We revisit the relation between H 2 maser detection rate and nuclear obscuration for a sample of 1 14 Seyfert galaxies, drawn from 
the CfA, 12^jm and IRAS F25/F60 catalogs. These sources have mid-infrared spectra from the Spitzer Space Telescope and they 
are searched for X-ray and [O III] 5007 A fluxes from the literature. We use the strength of the [O IV] 25.9yum emission line as tracer 
for the intrinsic AGN strength. After normalization by [OIV] the observed X-ray flux provides information about X-ray absorption. 
The distribution of X-ray / [O IV] flux ratios is significantly different for masers and non-masers: The maser detected Seyfert -2s (Sy 
1.8-2.0) populate a distinct X-ray / [O IV] range which is, on average, about a factor four lower than the range of Seyfert-2 non-masers 
and about a factor of ten lower than the range of Seyfert-ls (Sy 1.0-1.5). Non-masers are almost equally distributed over the entire 
X-ray / [O IV] range. This provides evidence that high nuclear obscuration plays a crucial role for the probability of maser detection. 
Furthermore, after normalization with [O IV], we find a similar but weaker trend for the distribution of the maser detection rate with 
the absorption of the 7 /jm dust continuum. This suggests that the obscuration of the 7 pm continuum occurs on larger spatial scales 
than that of the X-rays. Hence, in the AGN unified model, at moderate deviation from edge-on, the 7 /jm dust absorption may occur 
without proportionate X-ray absorption. The absorption of [O III] appears unrelated to maser detections. The failure to detect masers 
in obscured AGN is most likely due to insufficient observational sensitivity. 
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1. Introduction 

H2O megamaser galaxies represent an extreme subclass of ac- 
tive galactic nuclei (AG N) wit h st rong water maser em ission at 
22 GHz (reviews bv lLoll2005l and iHenkel et al.ll2005l) . In those 
cases where the emission arises from a molecular disk and can 
be resolved spatially using Very Long Baseline Interferometry, 
the central black hole (BH) mass and the distance to the galaxy 
can be determined (e .g. for NGC4258, lGreenhilletal.Nl993l 
iHerrnstein et al.l fl999l> . Thus, finding megamasers (henceforth 
simply called masers) and understanding their properties is of 
great interest. 

From theoretical considerations, a large line-of-sight col- 
umn density of velocity coherent gas favors the detection of a 
maser. High velocity coherence of the maser emitting gas is re- 
quired, because energy and momentum conservation imply that 
the induced photon has t he same frequ ency and direction as the 
stimulating photon (e.g. lElitzurll20 02). While the emission of 
an individual maser spot is directional (i.e. beamed), a collec- 
tion of such spots statistically may be expected to radiate in 
all directions, but this has not been confirmed so far. The orig- 
inally discovered water maser emission from AGN comes from 
(presumably edge-on) disks, and the resolved emission in most 
sources traces accretion disks while a few cases are star forma- 
tion masers. However, two sources, Circinus and NGC 3079, 
show in addition also off-disk jet masers that seem to trace out- 
flows. Potentially the se outflow masers are actually torus clouds 
dNenkova et al.ll2008l) . 

In the AGN unified model, an optically thick obscuring dust 
torus is envisioned to encircle the accretion disk and ty pe- 1 AGN 
are seen pole-on while type-2 AGN are seen edge-on (lAntonuccil 



1991 . Masers are almost exclusively found in AGN of Seyfert-2 
or LINER type, consistent with the picture t hat masers are pref- 
erentially beamed in the plane of the torus dBraatz et al.l [l 997. 
2004 IHenkel et alf2 005). But not all type-2 AGN are masers. 

From the conceptional viewpoint, it should be noted that the 
22GHz radio-frequency maser emission itself is believed to be 
largely unaffected by absorption; but a high X-ray, optical or 
mid-infrared obscuration may signpost a high likelihood that the 
masing disk is seen edge-on, hence favoring a maser detection. 

Type-2 AGN that host masers show a prevalence (> 80%) 
of high X-ray obscuring columns (7V H > 1 23 cm 2 ) and about 
half are Co mpton thick (N H > 1 24 cm 2 ) dBraatz et all 1 19971 
Zhang et al.ll20Q6l iG reenhill et al. 2008). However, as pointed 
out by IZhang et all d2006l) . among type-2 AGN the average X- 
ray derived column densities of masers and non-masers 1 are in- 
distinguishable. One explanation for this unexpected result could 
be that X-ray scattering in clumpy media dilutes the true line-of- 
sight column density, and thus prevents us from deriving unbi- 
ased orientation information. Therefore it is vital to include also 
information from other than X-ray wavelengths, to reveal the po- 
tential influence of nuclear obscuration on the maser detection 
and non-detect ion, respectiv ely. 

Recently, Zhan g etakl d2010h analyzed the Ka iron 
line equivalent wi dth EW(Kor), following the strategy of 
Bassani et alj d!999l) . and compared it with two optical thick- 
ness parameters, the infrared 6-400 pm luminosity Lir derived 
from IRAS 1 2-1 00 pm photometry and the [O III] 5007 A emis- 
sion line luminosity L[oni]- Both parameters were adopted to 



We denote as non-masers those AGN that have been observed at 
22 GHz, but for which no megamaser was detected. 
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be isotropic tracers for the intrinsic AGN strength. While the 
EW(Ka) distributions of 19 masers and 34 non-masers cover 
the same broad range (100 - 3000 eV), the median EW(Ka) of 
masers is about a factor 4 higher than that of the non-masers, 
indicating that the X-ray continuum of masers is more absorbed 
than that of non-masers. 

However it is still a matter of debate, whether Lir and 
L[oiii] are indeed isotropic tracers of the intrinsic AGN lu- 
minosity. While rOHIl has often been used as isotropic 



AGN tracer (iMulchaev et al.lfl 994; 



Bassani et al. 1999; Heckman et al 



Alonso-Herrero 



2005; Panessa 



et all 
et all 



1997 



2006; 



Lamastra et al. 120091) . the discovery of polarized [O III] emission 



in some type-2 AGN ddi Serego Alighieri et al.lll997l) cautions 
that a substantial fraction of [O III] can be shielded by the torus. 
Further studies, using MIR emission lines like [O IV] or [Ne V] 
as orientation independent tracer of the AGN power, provide ev- 
idence by means of the [OIII]/[OIV] ratio, that [OIII] suffers 
orientatio n dependent extinction, up to a factor of 1 in individ- 
ual ca ses (lHaas et al. 1 120051 Melende z et alj|2008al iBaum et all 

l2oToh . 

This is qualitatively consistent with results obtained us- 
ing the (extinction corrected) 2- lOkeV X-ray lumin osity Lx 
as intrinsic AGN power measure; iNetzer et al.l (120061) find that 
L[Oiii]/Lx of type-2 AGN is, on average, about a factor two 
lower than that of type-1 AGN. From a conceptional viewpoint, 
even in the face-on Syl case, the back-sided cone of the NLR 
lies - at least partly - behind an absorbing layer (e.g. the dust 
torus). Therefore it is highly questionable how far [OIII] can 
serve as an isotropic AGN tracer. The extinction correction via 
B aimer decrement (H a / Hp — 3) remains highly uncertain, since 
it is dependent on the geometry of the emitting and obscuring 
reg ions. 

iHes et alj d!996l) ; iBakerl d 19971) caution against the use of 
[O III] as a measure of the intrinsic NLR emission and suggest to 
use [Oil] 3727 instead. Observations of radio-loud AGN, where 
the orientation can be inferred from radio morphology, show that 
[Oil] is largely orientation independent (IHes et al.ll 19961: [Baker 
|1997|) . On the other hand, because of its low ionization potential, 
[OH] can also be dominated by star formation in the host (e.g. 
lHoll2 005). Then, the d e cline of L r p tti /L x with increasing L x , as 
found bv lNetzer et al.l (120061) . could be naturally explained by a 
decline of host/ AGN with increasing AGN Lx- 

Likewise the mid-infrared (A < 40/um) part of Lir is orienta- 
tion dependent (e.g. Fig. 16 in Buchan an et al.1 12006). while the 
far-infrared (A > 40/im) emission of Seyfert galaxies and low- 
luminosity quasars actually is dominated by star forming con- 
tribut i ons rather than by the A GN itself (e.g. iMaiolino & Riekd 
1 1 9951 ISchweitzer et~all r2006). Thus, a careful re-investigation 
using more suited isotropic AGN tracers would be desirable. 

Here, we revisit the connection between maser detection 
rate and nuclear obscuration using the strength of the [OIV] 
25.9/mi emission line (for short [OIV]) as tracer for the in- 
trinsic AGN strength. [OIV] has been found to be largely 
unaffected by obscuration (e.g.. iGenzel et al]|1998l [Haas et al. 
120051 Melende reTal.l2008a).lBaum et al.l2010h. We com bine the 
strategies of iBassani et al.l (1 19991) and lMelendez et al.1 d2008a). 
The observed X-ray (2-10 keV) flux normalized by [OIV] 
should provide information about X-ray absorption, even in 
the case of X-ray scattering caused by a complex geometry or 
for Compton thick cases. We compare the distribution of X- 
ray / [OIV] for masers and non-masers. In addition, after nor- 
malization with [OIV], we inspect the relation between maser 
detection rate and absorption of the 7 /jm dust continuum emit- 
ted from the nuclear torus, as well as maser detection and the 



absorption of the [O III] 5007 A emission of the central part of 
the narrow-line-region (NLR). 

The distances from which we derived the luminosities are 
taken from the NED database. The cosmology is based on H D = 
73 km s- 1 Mpc- 1 , Q A = 0.73 and Q,„ = 0.27. 

2. Data 

2. 1 . The parent sample 

At first glance, one could take all known masers and non-masers 
from the literature and compare their properties, for instance 
Lx/L[oiv]- But in order to determine nuclear obscuration, one 
needs to know also the range of Lx / L[o ivj for unobscured (pref- 
erentially Syl) sources, which should comprise a complete sam- 
ple free from any selection bias. However, the list of Seyferts, for 
which a maser search has been performed, did not follow clear 
selection criteria. Even worse, most maser searches have been 
performed on Sy2s, but only a small number on Syls. Because 
incomplete sample selection may influence the results, we here 
decided to start with complete Seyfert catalogs having well de- 
fined selection criteria. In order to increase the sample size, we 
created a master sample from the following three catalogs, con- 
sisting of a total of 163 sources. 

• The magnit ude limited complete s ample of the CfA Redshift 
Survey bv iHuchra & Burgl (fl992). which was supplied with 
updated Seyfert-type information from the NED database. 

• The 1 2//m Active Galax y Sample bv ISpinoglio & Malkanl 
(I1989I) . complemented bv lRush & Malkanl (1 1993h ~ 

• The IRAS F25/F 60 flux-ratio selected sample by 
Ide Griip et ail (Il992l) . as refined bv lSchmitt et al.l (120031) . 

Table 1 documents how the 163 sources distribute over the three 
catalogs, and how these catalogs match or complement each 
other. In general, we will present the results for the combined 
sample, but - where necessary - also for the catalogs individu- 
ally (Tab.©. 

Further below (Section 13.4) we will discuss potential dif- 
ferences between the three samples and our combined sample 
and all other known masers outside of it. The Spitzer data 
archive contains IRS spectra (at ~ 26/mi) for 126 of the 163 
sources classified as Seyferts according to the NED. This data is 
listed in Tab. [3] It covers the complete CfA-sample of 54 Seyfert 
Galaxies. It includes 107 of 1 18 Seyferts (two Blazars included 
as Syl) from the 12/im selected sample. For the IRAS sample 
we found useful IRS spectra for 34 of 60 sources. 

2.2. Maser information 

The parent sample of 126 sources with Spitzer spectra was 
searched for known maser-detections and non-d etections. For 
this p urpose we used the lists as compiled by iBennert et al.l 
(2009) and on the website of the Hubble Constant Maser 
Experiment (HoME) 2 . 

This search results in 18 masers (3 Syls, 15 Sy2s), 96 non- 
masers (36 Syls, 60 Sy2s) and 12 sources (10 Syls, 2 Sy2s) for 



http://www.cfa. harvai'd.edu/~lincoln/demo/HoME/surveys/survey.html 
compiled from Kondratko et al. 1 2006b a); Haeiwara et al. 12003); Zhang et al. 
(2006); [ Braatz et alj 120041) ; lHagiwara et al l d2002h; IGreenhill et all 42002L 
M997I); Isato et alj 120051); iBraatzetalJ <1996D; IGreenhill et alj 419951); 
iNakai et all f l 995); Hen kel et 311 )1998. 2005); B raatz et alJ<2003h:lHenkel et alj 
<1984D;lHaschick &Baari )198a);IClaussen & Lol )1986l):lHenkel et alj 119861); 
iBecker et alJ)1993l);IGreenhill et al.Ul990t) ; lBraatz & Gugliuccil )2008l) : IHeiikel 
(2008); Braatz (2008) 
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which no maser search has been performed so far (henceforth 
called maser-unknown). The results are listed in Tab. [3] Col. 2. 

2.3. [O IV] 25.89 /jm Line and 7 jim Continuum Flux 

Our analysis is based on public archival IRS spectra of Seyfert 
galaxies. We used the post-basic-calibration data (PBCD), as re- 
duced by the Spitzer Science Center's (SSC) pipeline. This in- 
cluded droop-, stray-light-, cross-talk- and saturation correction, 
dark subtraction, flatfielding and coaddition. 

If possible, the IRS high resolution spectra with R ~ 600 
have been chosen, to avoid contamination of [OIV]25.89// m with 
the neighboring [Fell] 25.99 Mm emission line. If high-resolution 
spectra were not available, the low-resolution spectra were used, 
including a background subtraction, which was also performed 
by the SSC pipeline. 

For the high resolution data, collected with the shorter (4.7 x 
11.3 and 1 1.1 x 22.3 arcsec) slits, separate background observa- 
tions had to be chosen to evaluate the background contribution. 
This was performed in iRamolla (120091) . by comparing the back- 
ground with the source fluxes at the presumably weakest part of 
the source spectrum between 9 and 10/mi rest frame; with the 
result that the background contribution is negligible in compar- 
ison with the conservatively assumed flux calibration errors of 
15%. The resulting errors are calculated from an assumed 15% 
flux calibration error and the error of the line fitting routine. 

The [O IV] flux has been extracted by fitting a simple spec- 
tral model in a wavelength- window of ~ 0.3 //m around the 
[O IV] line. This model consists of a linear base, convolved with 
Gaussian profiles that also include the neighboring [Fe II] line. 
No [OIV] aperture corrections had to be applied, because for 
both, high- and low-resolution data, the slit apertures cover an 
area larger than the expected size of the NLR, as estimated from 
the relationship 3 of [OIII] luminosity to NLR size, found by 
iBennert et all d2002l) . 

We calculated the 7.6 yum (henceforth for short 7 fim) con- 
tinuum flux from the background-subtracted IRS low -resolution 
spectra. We used a modified version of the PAHFIT code by 
Smith et aTJ (120071) which estimate d the continuum in the 5 to 
11.8 jum branch. As suggested by Smith & Drainel (120081) . we 
did not correct the continuum fit for the silicate feature at 9 .7 fim. 
The 7 jum continuum flux is then calculated from the PAH flux 
and the equivalent-widt h of the features at 7.4 /mi, 7.6 fim and 
7.8 yum (IRamol la 2009). The uncertainties of the 7 /um contin- 
uum are conservatively estimated to be smaller than 30%, which 
is sufficient for our purpose. In a few cases the AGN contri- 
butio n may be conta minated by nuclear (<3.7") star formation 
(e.g. lDeo et ail 2009). We checked that the effect on our statis- 
tical analysis is negligible by comparison with high resolution 
ground-based MIR observations. 

The [O IV] and 7 //m fluxes are listed in Tab. [3] The val- 
ues are consistent with those derived by others ([B uchanan et all 
l2006llDeo et al]l2007llDeo et al.ll2009llT~ommasin et alJl2010h . 

2.4. X-rays + [OIII] 5007 A line from the literature 

The 2-10keV hard X-Ray data have been obtained by several 
observers using ASCA, Beppo SAX, Chandra and XMM. We 
collected the data from the NED; in case of multiple entries we 
chose the latest detection. 

We have collected [O III] 5007 A emission line fluxes from 
various literature sources, as listed in Tab. [3] Because of the 

3 tog(R NLR ) = (0.52 ± 0.06) x log(L [0m] ) - (18.5 ± 2.6) 



large uncertainties, we did neither apply any extinction correc- 
tion nor any aperture correction for the [O III] fluxes. Such aper- 
ture corrections would affect a few very nearby sources, but most 
sources are sufficiently distant so that in the statistical analysis 
any bias is small. 

2.5. Additional maser sources 

On the one hand, our combined sample is drawn from the CfA, 
12yi/m and IRAS F25/F60 catalogs, containing 15 Sy2 maser 
sources with Spitzer spectra. On the other hand a total of 52 
masing Sy2 are known so far (IBennert et al .1120091 IZhang etlril 
2010), although drawn from different AGN catalogs using inho- 
mogeneous criteria. 

In order to compare the 15 Sy2 maser of our combined sam- 
ple with the remaining 37, we also analyzed available Spitzer 
spectra and gathered further [O III] and X-ray fluxes for them 
from the literature in the same manner, as we did on our com- 
bined sample. This results in an "off-sample" list of 37 Sy2 
masers that is appended to Tab. [3] 

3. Results and discussion 

While our combined sample of 126 sources contains 1 14 maser 
and non-maser sources (12 maser-unknown excluded), not all 
of them have data in all observables considered here (X-rays, 
F7/im> [OIII]). Therefore, we compare the maximum possible 
subsets for pairs of observables, and discuss the implications in 
the framework of the AGN unified model. Therein we consider 
as components the accretion disk, supposed to house the maser, 
the dust torus, the bi-conical NLR and the host galaxy. We here 
denote by Syls the subtypes between Seyfert 1.0 and 1.5, and 
by Sy2s those between Seyfert 1.8 and 2.0. All Figures con- 
tain an combined error bar in the lower right corner that is aver- 
aged from all relative errors in this measure. Since the literature 
sources did not uniformly presented errors, we do not perform 
this step for the X-ray, [OIII] and H2O luminosities. 

3. 1 . Nuclear X-ray obscuration 

Figure \T\ shows the observed hard X-ray versus [O IV] line lu- 
minosities, and Figure |2]the Lx/L[oivj histogram. The features 
are: 

• On average, Syls present an about 10 times higher X- 
ray/ [OIV] ratio than Sy2s. 

• Sy2 non-masers are evenly distributed over the entire range 
occupied by Sy2 masers and by Syls (Fig.[TJ. 

• Sy2 masers are almost disjoint from Syls. At a given [OIV] 
luminosity, Sy2 masers have on average about a factor 
4 lower X-ray luminosity than Sy2 non-masers (Fig. |2}. 
Likewise, the few Syl masers have a lower X-ray luminosity 
than Syl non-masers. 4 

• 6 out of 12 sources without masing information, but Lx and 
L[oiv] available (see Tab. [3}, show the same trends as the 
Syls and Sy2s with masing information (see Tab. 0. They 
are not plotted, to avoid overcrowding of Figs. [T]and|2]with 
too many different symbols. 



4 The Syl masers are NGC4051, NGC4151 and UCG5101. Note 
that both NGC4051 and NGC4151 have relatively low maser H 2 lu- 
minosity, and UCG 5 101 is an ultra-luminous infrared galaxy so that the 
maser luminosity could arise from starburst regions rather than from the 
AGN accretion disk. 
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Fig. 1. Observed 2-10 keV X-ray versus [OIV] line luminosity 
Blue stars represent Syls (Sy 1.0-1.5), red circles Sy2s (Sy 1.8- 
2.0). Filled symbols are masers, open symbols are non-masers. 
The dotted lines mark fixed Lx/Lpiv] ratios of 1000; 100; 10; 
1; 0.1 (from top to bottom). The error-bar in the lower right 
corner is the average relative error of all [O IV] measurements. 
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Fig. 2. Histogram of the Lx/L[oivj ratio of the data points 
shown in Fig.Q] The red dashed line represents Sy2 non-masers, 
whereas the blue thick line represents the Syl non-masers. The 
maser-detections are represented by the dashed surfaces - blue 
and vertically dashed for Syl, red and horizontally dashed for 
Sy2. The one upper limit is excluded. 



We assume that the X-ray deficit, i.e. the decrease of Lx/L[oiv], 
is caused by obscuration, probably in the molecular dust torus. 
Then, the Figures Q] and [2] clearly demonstrate that masers are 
found almost exclusively in Sy2s with heavy nuclear obscura- 
tion, while non-maser Sy2s exhibit a broad range of X-ray ab- 
sorption. A two-sided Kolmogorov-Smirnov (KS) test, applied 
on the Lx/L[oiv] distribution, results in a probability of 18% 
that the Sy2 masers and non-masers are drawn from the same 
parent population. Our results agree with those of Greenhill et al. 
(120081) and lZhang etal] d2010h who find that about 60% of the 
masers are Compton-thick. 

Assuming that Syls are almost unobscured, the obscured 
sources populate the Lx/L[oiv] range below 10 in Fig Q] 
Thus, masers populate almost completely the range of obscured 
sources. Surprisingly this range also contains numerous non- 
masers. In order to better understand why in such absorbed 
sources the maser search failed, we consider the influence of 
observed brightness. Fig. |3]shows the X-ray and [O IV] flux dis- 
tribution (instead of the luminosity distribution). In fact, the Sy2 
masers and non-masers show a flux-dependence in their [O IV] 
distribution. Sources with low [OIV] flux are more frequently 
classified as non-masers (2 Sy2 maser and 12 Sy2 non-maser 
at L [0 iv] < 1CT 16 WirT 2 ), while sources with high [OIV] flux 
are more frequently classified as masers (6 Sy2 maser and 5 Sy2 
non-maser at L [0 iv] > 10~ 15 WirT 2 ). Fig. |4] displays the [OIV] 
fluxes of all Sy2s lying below the dividing line between obscured 
and unobscured sources (Lx : L[oiv] ~ 10). Among this subset 
of obscured Sy2s, the frequency of non-masers rises constantly 
towards lower [OIV] fluxes, in contrast to the distribution of 
masers. This incidence is consistent with an observational bias 
against the maser-detection for faint AGN. This implies a rela- 
tion between [O IV] flux and H2O flux which is indeed observed 
in Fig. [9] and discussed in Section l3~31 

Because the detection of maser emission appears to be biased 
against sources with low flux, we conclude that among obscured 
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Fig. 3. Observed X-ray versus [OIV] line flux. Symbols and 
colors are as in Fig. Q] The dotted lines mark fixed flux ratios of 
1000; 100; 10; 1; 0.1 (from top to bottom). 

The error-bar in the lower right corner is the average relative error of all 
[OIV] measurements. 



sources the true fraction of masers is higher than indicated by 
Fig.El 



3.2. Extended obscuration of the dust torus and the NLR 

Figure [5] shows a histogram of the MIR 7 /im continuum to 
[O IV] line ratio. The striking features of this diagram are: 
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Fig. 4. [O IV] flux histogram of Sy2 masers for absorbed sources 
from Fig.[3]wifh Lx / L[oivj < 10. Masers are represented by the 
shaded area, non-masers by the thick dashed histogram. 
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Fig. 5. Histogram of the 7 /mi continuum to [OIV] line ratio. 
Legend as in Fig. [2] All 21 upper limits from Tab. [3] are ex- 
cluded. 



• Sy2s populate about the same total range as Syls, but show 
a prevalence for lower 7 /mi / [OIV] values, i.e. a 7 /mi 
continuum deficit. On average, the ratio 7 /mi / [OIV] of 
Sy2s is about a factor of 3 lower than that of Syls. This is 
consistent with the results obtained via radio normalization 
(7 jum / 8GHz) by Buchanan et al. (2006, their Fig. 16). 

• Among Sy2s the 7 /mi / [O IV] ratio of masers is, on average, 
about a factor of 2 lower than that of non-masers. A KS test 
results in a probability of 3.7% that the Sy2 masers and non- 
masers are drawn from the same parent distribution. Flux 
considerations similar to those for Lx / L[oiv] suggest that 
the true 7 /mi / [O IV] separation of masers and non-masers 
will be even more pronounced once the observational bias 
against the detection of low flux masers is taken into account. 

We assume that the deficit of the 7/mi continuum in Sy2s is 
mainly caused by absorption of the torus dust emission. This ab- 
sorption has to take place somewhere between the emitting re- 
gion and the observer, hence probably in the "halo" of the torus, 
i.e. in the outer part of the torus itself or in the host galaxy. It 
is possible that the scale height of this MIR-absorbing halo, i.e. 
the projected distance of absorbing material from the line-of- 
sight to the nuclear accretion disk, is (much) larger than the scale 
height of the torus itself. This is consistent with the results from 
a Spitzer study of CfA Seyferts (iDeo et al .1120071) . where sources 
with high 10 //m silicate absorption show a preference for large 
host-galaxy inclinations and irregularities (merger events or in- 
teractions), both of which lead to absorption through the host. 

In order to provide further clues on the extent of the MIR- 
absorbing material, we consider the [OIII] 5007 A versus 
[O IV] line luminosity as shown in Figures [6] and [7] The fea- 
tures of the [O III] / [O IV] distribution are similar to those of 
7 fim I [OIV]. Most Sy2s populate the same range as Syls, a few 
Sy2s show a [O III] deficit, i.e. on average about a factor 3 lower 
[OIII] / r OIV] ratios when compared to Syls, consistent with 
results by Bau m et aD(l2010l) on the 12 /mi sample. The distribu- 
tion pattern of masers and non-masers appears to be statistically 
indistinguishable. A KS-test results in a probability of 61% that 
the Sy2 masers and non-masers are drawn from the same par- 
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Fig. 6. Distribution of [OIII] 5007 A versus [OIV] 25.9/mi lu- 
minosity. Symbols are as in Figure Q] The dotted lines mark 
fixed-ratios 10; 1; 0.1; 0.01 (from top to bottom). The error- 
bar in the lower right corner is the average relative error of all 
[OIV] measurements. 

ent population. However, the Sy2 subsample in Figures [6] and [7] 
shows a distinct tail towards lower ratios, potentially caused by 
absorption. 

The [OIV] 25.9/mi line is ~50 times less affected 
by extinction than the optical [OIII] 5007 A line. A low 
[O III] / [ O IV] ratio argues i n favor of large obscuration, as ex- 
plained in lHaas et alj (12005). Another explanation of deviating 
[O III] / [O IV] ratios could be different radiation fields in the 
NLR. Because the [OIV] 25.9 /mi line needs a higher ioniza- 
tion potential than the optical [OIII] 5007 A line, AGN with a 
hard radiation field are expected to show a low ratio. The Sy2s 
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Fig. 7. Histogram of the [O III] to [O IV] ratio. Legend as in Fig 
|2] Five upper limits are excluded. 



with low [OIII] / [OIV] would then be those AGN with hard 
radiation fields. But this is not consistent with other spectro- 
scopic MIR tracers like the [Nell] 12.8 fim to [OIV] flux ratio 
dMelendez et al. 2008b). Thus we conclude that in Figures[6]and 
Qthe NLR of Sy2s with [OIII] deficit is considerably obscured. 

While the [OIII] obscuration may occur mainly in the in- 
nermost part of the NLR, the large extent (>1 kpc) of the NLR 
suggests that the absorption is not confined to the region en- 
circled by the dust torus. Rather the sky-projected distribution 
of the absorbing material might reach further out to a consider- 
able distance (several hundred parsec) from the line-of-sight to 
the nuclear accretion disk. The presence of moderately extended 
[O III] absorption, as well as the similarity of the 7 fim / [O IV] 
and [O III] / [O IV] diagrams supports the picture that also the 
MIR absorption takes place in a moderately extended layer, i.e. 
the torus halo mentioned above. Although both observables 7 fim 
and [O III] appear to be affected by absorption in a similar fash- 
ion, we note that the distribution of 7 fim / [O III] spans a large 
range (2-700). This is not surprising in view of the diversity of 
the orientation-dependent appearance of the involved emitters 
and absorbers even for a simple AGN model. 

3.3. Combined picture 

Why do masers and non-masers show so different distributions 
in X-ray / [O IV], while their distribution in 7 fim / [O IV] looks 
more similar? 

Because masers need a large line-of-sight column density of 
velocity coherent gas, they are expected to be predominantly de- 
tected in edge-on accretion disks. Thus, the maser detection or 
non-detection can tell us about the disk orientation with respect 
to the line-of-sight. In order to constrain the implications in the 
framework of the AGN unified model, we consider two extreme 
cases: 

1. For a disk seen edge-on, obviously the maser is most eas- 
ily detected and the molecular torus is seen more or less 
edge-on, too. In this case the nuclear accretion disk (and its 
corona) is shielded by the torus, so that the X-rays are heav- 
ily obscured. If additional extended material, able to ob- 



scure the MIR emission, does not lie in the torus plane, the 
7 fim I [O IV] ratio is decreased. 
2. For a non-maser, both, disk and torus appear to be suffi- 
ciently tilted away from edge-on, so that the nuclear X-ray 
absorption is relatively low. In addition, our diagrams indi- 
cate the existence of non-masers, where the torus plane is 
seen edge-on, but the disk could be tilted out of this plane 
due to locally different angular momentum. In this case of 
a non-maser, the edge-on torus causes a high obscuration of 
the X-ray nucleus as well. On the other hand, irrespective 
of the disk and torus orientation, the MIR continuum can be 
absorbed or not depending on the line-of-sight through the 
extended host. 

From these two extreme cases we see: While the requirement 
for heavy absorption of the nuclear X-rays is that the line-of- 
sight has to hit a rather compact area with very high column 
density, the area of the torus emission and even more the area of 
the (bright) NLR emission is orders of magnitude larger, so that 
the absorber must cover a larger area, too. If the absorption of 
X-rays and 7 fim occurs on different spatial scales, the strength 
of the obscuration in each wavelength range may be sensitive to 
small differences in the aspect angle. Furthermore, because the 
MIR-absorbing material is located farther away from the line-of- 
sight to the nucleus, it is less reliable to predict whether a maser 
will be detected. 



3.4. Comparison of the three samples with other known 
masers 

Our combined sample was compiled from three complete sam- 
ples with good coverage in the Spitzer archive. Table Q] shows 
the overlaps between the samples. Note that each sample is in- 
complete due to the limited availability of data in the Spitzer 
IRS archive and of X-Ray and [O III] measurements in the lit- 
erature (see Tab. [2). Moreover, maser surveys were not per- 
formed with homogeneous properties (sensitivity, velocity cov- 
erage) and were not carried out for all sources of our combined 
sample. 

The three samples were based on different selection criteria: 
Optically selected Seyferts in the CfA sample and IRAS selected 
sources in the 12yum and F25/F60 sample. Thus, it is possible that 
they suffer from different biases with respect to potential maser 
detection. The fraction of Sy2 masers to non-masers increases 
from 1 1 A (4/16) in the CfA sample, to ~ 1/3 (1 1/29) in the 12jum 
sample and to > 1 (8/6) in the IRAS F25/F60 sample. This is 
consistent with the well known fact that the mid- and far-infrared 
wavelengths select more obscured AGN than the optical bands. 

However, the range of luminosity ratios (Lx / L[orv], 
Lv^m / L[oivj, L[oni] / L[oiv]X listed in Tab. [2] are similar for 
all three samples. This indicates that also among optical selected 
masing sources, some can be obscured with a level, similar to 
that of infrared selected sources 5 . To summarize, the result of 
all three samples (CfA, \2fim and IRAS F25/F60) are similar in 
that they point consistently to a prevalence of maser detections 
in Sy2s with high X-ray obscuration and one may expect this 
holds also for Seyfert galaxies in general. 

Are our s elected Sy2 masers representative for all 52 known 
Sy2 masers dZhang etalJ 12010b ? To address this question, we 
compare our in-sample Sy2 masers with all remaining 37 off- 
sample Sy2 masers. 



5 Even for infrared selected AGN, optical criteria influence the sam- 
ple, because the Seyfert identification is done by optical spectroscopy. 
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In Figure [8] we show a comparison of X-ray and [OIV] lu- 
minosities between in- and off-sample masers. The compari- 
son refers to those masers with X-ray and [OIV] fluxes avail- 
able, i.e. 12 off-sample Sy2s 13 in-sample Sy2s and 3 in-sample 
Sy Is. Both, in-sample and off-sample roughly populate the same 
Lx / L[oiV] range. But the Lx / L[oiv] ratio is, on average, about a 
factor 2 higher for the off-sample than for the in-sample masers. 
This indicates that the off-sample Sy2 masers may be less ab- 
sorbed than the in-sample ones. Compared with the in-sample 
Sy2 non-masers (omitted in Fig [8] see Fig. [TJ, however, the off- 
sample masers show, on average, still about a factor 2 lower 
Lx / L[oiv], hence considerably high obscuration. 

Some off-sample masers show spurious flux ratios that im- 
ply no obscuration (i.e. Lx / Lpiv] > 10 for 4 objects). 
Among them, we find two nearby extended sources, NGC4258 
and NGC4945, in which t he X-ray emiss ion has been associ- 
ated w ith star formation by Risaliy (120021) and IStrickland et al.l 
J200l . 

To summarize, the large overlap and the lack of significant 
differences between in- and off-sample Sy2 masers indicates that 
the results for our combined sample's Sy2 masers can be ex- 
tended to all known Sy2 masers. 

We note that the inhomogeneous selection of all off-sample 
masers and non-masers precludes to derive a meaningful com- 
parison of maser to non-maser statistics with our in-sample data. 
A KS test shows a probability of 63% that both subsets, off- 
and in-sample masers, are drawn from the same parent distri- 
bution and the above mentioned difference in Lx / Ljoivj is 
only by chance. But similarly a KS test between the in-sample 
non-masers and the off-sample masers yields also 60% prob- 
ability that they are drawn from the same parent distribution. 
But yet, comparing the Lx / L[oiv] ratio between in-sample Syl 
non-masers and the off-sample Sy2 masers, shows a probability 
of 0.16%o to be drawn from the same parent distribution. This 
shows that the off-sample Sy2 masers are still significantly dif- 
ferent from the unobscured Syl non-maser. 

3.5. Maser and AGN luminosity 

A search for H2O masers in 274 high-redshift (0.3 < z < 
0.8) SPSS type-2 A GN half of which being type-2 quasars 
(iBennert et alJl2009h found on ly one maser (SDSSJ0804+3607, 
Barvainis & Antonucci 2005). The high rate of non-detections 
in these luminous AGN could be due to limited observational 
sensitivity or to intrinsic differences between low- and high- 
luminosity AGN. Such differences could be, for instance, that 
in a high-luminosity AGN the accretion disk becomes hotter so 
that the density required for maser emission falls below a critical 
limit. If this is frequently the case, one would expect a relative 
decline of H2O maser luminosity with increasing AGN luminos- 
ity. On the other hand, the SDSS H2O maser survey was rela- 
tively shallow, because one was interested to find masers which 
are sufficiently bright for spatially resolved follow-up VLBI ob- 
servations. 

Here, we consider how far the Seyfert sample can help to dis- 
tinguish between these two possibilities (i.e. by looking whether 
or not Lh20 / L[oivj declines with increasing Ljoiv])- A remark- 
able feature of Fig.[T]is that maser-detections and non-detections 
are quite evenly distributed along the whole [O IV] luminosity 
range covering about 4 orders of magnitude. Thus our data do 
not indicate a trend that the frequency of non-masers rises with 
luminosity. We have also seen that the available maser observa- 
tions of the Seyferts are biased against maser detection in faint 
(and distant) AGN (Fig. g}. 
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Fig. 8. Observed 2-10 keV X-rays plotted against [OIV] line 
luminosity. Syl masers are represented by blue stars and Sy2 
masers by red dots. Off-sample Sy2 masers are marked with 
black dots. The dotted lines mark fixed flux ratios of 1000; 100; 
10; 1; 0.1 (from top to bottom). The error-bar in the lower right 
corner is the average relative error of all [O IV] measurements. 
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Fig. 9. Maser H2O versus [O IV] luminosity. Symbols and col- 
ors are as in Fig. [8] The dotted lines mark fixed-ratios 10~ 4 ; 
10" 5 ; 10" 6 ; 10" 7 (from top to bottom). For comparison, the ex- 
pected position of the z = 0.66 maser SDSSJ0804+3607 is 
marked with a black dot and a gray bar covering the range 
0.1xL[oih] < Ljoiv] < 100 xL[o in] assumed from Fig. |6] 



Figure [9] shows the (isotropic) maser luminosity versus 
the AGN luminosity as traced by [OIV]. In addition to the 
Seyfert sample we have plotted the expected position of 
SDSSJ0804+3607 at z = 0.66, the only QSO-2 maser detec- 
tion. Because this source was not observed with Spitzer IRS, 
we derived L[oivj from Lpm] using the range 0.1 x L[omj < 
L[oiV] < 100 x L[oni] as indicated in Fig. [6] which is also 
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valid for higher-luminosity AGN (Haas et al 1 120051) . At a given 
[OIV] luminosity the maser luminosity spreads over three or- 
ders of magnitude (Fig.|9j. One explanation for the large spread 
is that the maser emission is, in fact, not isotropic and hence 
the derived maser luminosity depends sensitively on the maser 
direction with respect to the line-of-sight. The Seyfert sample 
alone indicates only a marginal correlation in Fig. [9] with a 
Pearson correlation coefficient of 0.32 for all Sy2 masers of our 
combined sample that is not significant at the 5% level. Adding 
the off-sample Sy2 masers changes the coefficient to 0.46 which 
would then be significant, but the correlation could also be an 
artifact of distance in the luminosities. However, combined with 
the position of SDSSJ0804+3607 and the fact, that its assumed 
H2O / [O IV] ratio lies in the same range as for the lower lumi- 
nosity AGN, argues in favor of a physical connection between 
maser and AGN luminosity. 

The numerous non-masers among the SDSS QSO-2s 
dBennert et al.l l2009h have a [OIII] luminosity similar to that 
of SDSSJ0804+3607, hence are expected to populate a sim- 
ilar L[oiv] range in Fig. |9] The H2O maser upper limits 6 , 
found for these QSO-2s bv lBennert et al.l (120091) . lie even above 
J0804+3607. 

Thus the upper limits are not stringent enough to support a 
relative decline of maser luminosity with increasing AGN lumi- 
nosity. This, together with the sufficiently high Lh20 / L[oiv] 
ratio of J0804+3607, leads us to conclude that the main reason 
for the high rate of maser non-detections is insufficient observa- 
tional sensitivity, rather than basic differences between low-and 
high-luminosity AGN for hosting a maser. 

4. Conclusion 

In order to understand the connection between H2O maser de- 
tection rate and nuclear extinction we used the [OIV]25.9^ m 
line and the 7 fim continuum flux from Spitzer spectra of a 
well-selected sample of 114 Seyfert galaxies, from the CfA, 
12yum and IRAS F25/F60 catalogs, for which a maser search has 
been performed. These data were then compared to hard X-ray 
and [OIII] 5007 A fluxes from the literature. We analyzed the 
data in the framework of the orientation-dependent AGN unified 
scheme, yielding the following results: 

1 . Comparing hard X-rays to [O IV] flux, Sy2s exhibit, on aver- 
age, an about 10 times lower X-ray to [OIV] ratio than Syls. 
Masers prefer X-ray absorbed sources (i. e. low Lx / Ljoivj 
ratios). Sy2 masers present on average about 4 times less X- 
ray flux normalized by [OIV] than non-maser Sy2s. This is 
consistent with geometric alignment of both the X-ray ab- 
sorber and the the maser emitting region in the accretion 
disk. Non-masers do not show a preference for strongly ab- 
sorbed sources. However, our data indicate an observational 
bias against faint sources, in the sense that more sensitive 
maser observations might reveal more absorbed sources to 
house a maser. 

2. Regarding the 7 /mi to [OIV] flux ratio we find that most 
Sy2s spread along the same range as Syls. However there are 
sources with a significantly lower ratio, rendering the Sy2s 
on average about 3 times lower than Syls. These cases can 
be explained by an extended dusty absorber that is covering 
the 7 jt/m emitting torus region. Maser-detections also show 
a preference for 7 fim absorbed sources, but with less signifi- 
cance than in the X-ray to [OIV] comparison. This suggests 

e ipa = 0039 x j_ x ( Dj\ 2 



that the geometric alignment of the MIR absorber with the 
maser emitting disk is not as perfect as the supposed align- 
ment of the disk with the X-ray absorber. 

3. The [OIII] to [OIV] flux ratio also indicates the presence 
of extended obscuration in some Sy2s that blocks the op- 
tical emission from the NLR. Masers and non-masers are 
distributed very similarly in the [OIII] to [OIV] plot. The 
fraction of maser-detections is not significantly higher for 
sources with such extended absorption. This leads us to con- 
clude that the matter distribution for the [OIII] absorber 
is not essential for the prediction of a maser detection. 
Moreover, it is possible that a substantial fraction of the ab- 
sorption of AGN emission could occur in extended regions 
outside the torus that are not necessarily aligned with the 
torus or AGN sub-structure. 

4. The three samples, CfA, 12/mi and IRAS F25/F60 provide 
very similar results. The Sy2 maser to non-maser fraction 
increases from optical to infrared selection. The Sy2 masers 
of our combined sample have a similar range of L[oiv] and 
Lh20 as the known off-sample Sy2 masers. While, on aver- 
age, the off-sample masers are a factor two less obscured, as 
inferred by the Lx / L[oiv] ratio, they are still considerably 
obscured compared with Syls from our combined sample. 
Thus, the results obtained for our combined sample's Sy2 
masers may also hold for all remaining Sy2 maser sources 
that have no Spitzer or X-ray data available. 

5. After supplementing our combined sample with the remain- 
ing known Sy2 masers that were not included in it, the H2O 
maser luminosity appears to be correlated with the AGN lu- 
minosity as traced by [OIV], although with a large spread. 
We do not find evidence for physical differences between 
low-and high-luminosity AGN for housing a maser. 

The results demonstrate that heavy X-ray absorption is an indi- 
cator for high probability to detect a maser. The 7 /jm absorption 
can also be used to find maser candidates, but with lower proba- 
bility. 
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IRAS F25/F60 






60 






27 



Table 1. Documentation of how far the three samples (of 163 sources in total) match and complement each other. For each row the 
Table entries list the number of sources contained in the intersection of the sample denoted by the columns 2-7 with the CfA, 12/im 
and IRAS F25/F60 sample (Col. 1). E. g. the union of the 12/im and IRAS F25/F60 samples (Col. 5) has 45 sources in common 
with the CfA sample (Row 1). 
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Table 2. The average values and standard deviations of the logarithmic luminosit y ratios for each subset of Seyfert ga laxies. Column 
1: In descending order, the optically selected CfA sourc es (|Huchra & Burdl992h . the MIR selected 12/im sources dRush & M alkan 
[1991 . the IRAS F25/F60 flux-ratio selected sources (ISchmitt et al.l 120031) and the combined sample that is used in this work. 
Columns 2-4: Number of sources with luminosities available in [OIV] and one of the following: 2-10 keV X-rays (Col. 2) or 7/mi 
(Col. 3) or [OIII] (Col. 4). Each sample's row is subdivided into Syl maser, non-maser, maser-unknown and Sy2 maser, non- 
maser, maser-unknown. Column 5: Total Number of objects of the parent sample from Col. 1, e.g. 3 Syl maser in the CfA, 18 Syl 
non-maser in the CfA and so forth. 
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Additional off-sample data 
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Source 
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Table 3. Measured fluxes and literature values. Column 1: 

Source names with catalog marks: 

a: CfA sample 

b: 12/im sample 

c: IRAS F25/F60 sample 

Column 2: Isot ropic maser luminosities obtained from 
Ben nert et alj l2009i A question mark designates sources un- 
observed for masers and a dash represents a maser undetected 
source. 

Note that MRK93 8 and NGC1320 are listed as Maser in 
Zhan g et alj d2Q10t) but have no luminosity information avail- 
able. Therefore, they are marked with asterisks. Column 3: 
Seyfert type obtained from the NED or the literature with 
references give n in squ a re bra ckets. Column 4: [OIV] flux 
determined bv iRamollal d2009t) . Column 5: [OIII] flux from 
the literature. Column 6: 2-10 keV X-ray flux obtained from 
the literature. C olumn 7: 7 jjm continuum flux determined by 
Ramolla (2009). An analysis of the off-sample masers has also 
been performed. This data is appended in the table. 
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